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Abstract 


A  new  approach  for  design,  generation,  and  computerized  simulation  of  meshing  and  contact  of  face- 
milled,  formate  cut  spiral  bevel  gears  is  presented-  The  purpose  is  to  develop  a  low  noise,  stabilized  bearing 
contact  for  this  type  of  gear  drives.  The  approach  proposed  is  based  on  application  of  three  procedures  that 
permit  in  sequence,  to  provide  a  longitudinally  directed  bearing  contact,  a  predesigned  parabolic  function 
of  transmission  errors  and  limit  the  shift  of  bearing  contact  caused  by  errors  of  alignment.  The  theory 
developed  is  illustrated  with  an  example  of  design  and  computation. 
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Nomenclature 

Profile  angle  of  parabolic  blade  for  gear  head  cutter  at  mean  point  M 

Blade  angle  for  straight  blade  of  pinion  head  cutter 

Elastic  approach  of  contacting  surfaces 

Pitch  cone  angles  of  pinion  (/=7)  and  gear  (/=2),  respectively 

Machine  root  angles  for  pinion  (/=7)  and  gear  (/=2),  respectively 

Root  cone  angles  for  pinion  (i-I)  and  gear  (/=2),  respectively 

Shaft  angle 

Tangent  to  the  path  of  contact  on  the  pinion  (f=7)  and  gear  (/=2),  respectively 
Surface  parameters  of  the  pinion  and  gear  head  cutters,  respectively 
Surface  parameters  of  the  pinion  and  the  gear  fillet  parts  of  the  head  cutter 
Angle  formed  between  principal  directions  (ij,=  1,  2,  p,  g) 

Angles  of  rotation  of  pinion  (/=7)  and  gear  (/=2;  in  the  process  of  meshing 
Angle  of  rotation  of  the  cradle  in  the  process  for  generation  of  the  pinion 

Angle  of  rotation  of  the  pinion  in  the  process  for  generation 
Fillet  radii  for  the  pinion  and  gear 

Angular  velocity  of  the  pinion  (in  meshing  and  generation) 

Pinion  {i-I)  or  gear  (/=2)  tooth  surface 
Pinion  {k=p)  or  gear  {k-g)  generating  surface 
Pinion  and  gear  axial  displacements,  respectively 
Blank  offset  for  pinion 

Sliding  base  for  pinion 
Machine  center  to  back  for  pinion 
Machine  center  to  back  for  gear 
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Errors  of  the  offset  and  shaft  angle,  respectively 
Function  of  transmission  errors 
Parabola  coefficient  of  gear  head  cutter  blade 
Dedendum  of  pinion 

Mean  cone  distance  and  face  width 

Coefficients  of  modified  roll  for  pinion  generation,  (/=7,  2,  3) 

Unit  vectors  of  principal  directions  of  pinion  and  gear  tooth  surfaces,  respectively 
Unit  vectors  of  principal  directions  of  cutter  surfaces  for  pinion  and  gear,  respectively 
Normal  curvature  of  surface  E/ 

Relative  normal  curvature 

Principal  curvatures  of  pinion  and  gear  tooth  surfaces,  respectively 
Principal  curvatures  of  cutter  surfaces  for  pinion  and  gear,  respectively 
Cylindrical  coordinates  for  contact  point  A  on  gear  tooth  surface 
Second  derivative  of  transmission  function  (j>2i(l>i) 

gear  ratio 

Mean  contact  point 

Matrix  of  coordinate  transformation  from  system  Si  to  system  Sj 

Unit  normal  and  normal  to  the  surface  X  represented  in  coordinate  system  Sf 

Tooth  number  of  pinion  {i=l)  or  gear  (/=2) 

Installment  angle  for  pinion  head  cutter 

Radii  of  the  head-cutter  at  mean  point  for  the  pinion  and  gear 

Gear  cutter  radii 

Position  vector  in  system  5/  (/  =ly  2y  bj ,  b2  y  hy  k  mj  y  m2  y  py  g) 

Surface  parameters  of  the  pinion 

Coordinate  system  (/  =7,  2y  bj ,  b2  y  hy  k  mj  y  m2  y  p,  g) 

Radial  setting  of  the  head  cutter  of  the  pinion 

Gear  tooth  surface  parameter 

Velocity  of  contact  point  in  its  motion  over  surface 

Transfer  velocity  of  contact  point  in  its  motion  with  surface  E/ 

Relative  velocity  of  contact  point  (i  yj  =7,  2,  cj ,  p) 

Components  of  the  velocity  of  the  contact  point  in  its  motion  over  surface  E/ 

Center  distances  of  fillet  circular  arcs  of  pinion  and  gear,  respectively 
Length  of  major  axis  of  contact  ellipse 


NASA/CR— 2001-210894 


2 


1.  Introduction 


Design  and  generation  of  spiral  bevel  gears  is  a  subject  of  intensive  research  of  scientists,  designers  and 
manufacturers  [1  -  7],  High  quality  equipment  such  as  CNC  (computer  numerically  controlled)  machines 
and  tools  are  produced  by  the  Gleason  Works  (USA)  and  the  Klingelnber-Orlikon  (Germany  and 
Switzerland)  that  manufacture  these  gears.  However,  the  quality  of  spiral  bevel  gears  depends  substantially 
on  the  choice  of  machine-tool  settings  that  have  to  be  determined  for  the  generation  of  the  gears.  The 
machine-tool  settings  are  not  standardized  and  their  determination  is  the  subject  of  a  skilled  design.  The 
freedom  of  choosing  various  machine-tool  settings  provide  an  opportunity  for  generation  of  low-noise, 
stable  bearing  contact  spiral  bevel  gears. 

The  subject  of  this  paper  is  the  enhanced  design  of  face-milled  formate  cut  spiral  bevel  gears.  The  term 
“formate  cut”  means  that  the  gear  is  held  at  rest  where  it  is  cut  by  the  tool  (head-cutter).  Thus,  the  gear 
tooth  surface  is  a  copy  of  the  tool  surface  that  is  a  surface  of  revolution.  The  pinion  tooth  surface  is 
generated  by  a  head-cutter  while  the  tool  and  the  pinion  performed  related  rotations.  The  pinion  tooth 
surface  is  generated  as  the  envelope  to  the  family  of  the  tool  (head-cutter)  surfaces.  The  advantage  of 
formate-cut  method  is  the  high  productivity  of  gear  generation.  The  main  problem  of  the  design  of  formate- 
cut  spiral  bevel  gears  is  the  conjugation  of  tooth  surfaces  of  the  gear  and  the  pinion. 

The  contents  of  this  paper  cover  the  following  topics: 

(1)  Computational  procedures  of  proposed  design  that  provides  a  longitudinal  bearing  contact  and  a 
predesigned  parabolic  function  of  transmission  errors,  and  the  ability  to  investigate  the  influence  of 
misalignment. 

(2)  Computer  programs  and  algorithms  developed  for  synthesis  and  analysis  of  meshing  and  contact  of  gear 
drives. 

The  theory  developed  is  illustrated  with  an  example  of  computerized  design  of  a  formate  cut  gear  drive. 


2.  Basic  Ideas  of  Proposed  Approach 

Localization  of  Contact  [1,  8].  The  contact  of  pinion-gear  tooth  surfaces  is  localized  due  to  the  mismatch 
of  generating  surfaces  (surfaces  of  the  head-cutters)  that  is  achieved  by  application  of  different  dimensions 
of  the  gear-pinion  head-cutters. 

The  generated  pinion-gear  tooth  surfaces  are  in  point  contact  at  every  instant,  instead  of  in  line  contact.  The 
tooth  contact  under  the  load  is  spread  over  an  elliptical  area  and  the  bearing  contact  is  formed  as  a  set  of 
contact  elhpses.  The  magnitude  and  orientation  of  instantaneous  contact  ellipse  requires  information 
regarding  the  elastic  deformation  of  tooth  surfaces  and  the  principal  curvatures  and  the  directions  of  tooth 
surfaces  at  each  instantaneous  point  of  their  tangency. 

It  will  be  shown  later  that  the  path  of  contact  on  gear  tooth  surface  is  oriented  longitudinally.  Such  a  path  of 
contact  is  favorable  for  the  increase  of  contact  ratio  (the  average  number  of  teeth  that  are  in  contact 
simultaneously),  and  avoids  edge  contact  and  improves  the  conditions  of  transfer  of  meshing  where  one 
pair  of  teeth  is  changed  for  the  neighboring  one. 

The  TCA  (Tooth  Contact  Analysis)  computer  program  developed  by  the  authors  permits  simulation  of 
meshing  surfaces  and  the  shift  of  bearing  contact  caused  by  errors  of  alignment  [8,  9].  In  some  cases  of 
design,  it  may  become  necessary  to  adjust  the  bearing  contact  -  to  deviate  it  from  the  longitudinal  direction 
-  to  reduce  the  shift  of  bearing  contact  caused  by  misalignment. 
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Application  of  Parabolic  Profiles  of  Blades  of  Gear  Head- Cutter,  We  have  mentioned  above  that  the 
gear  tooth  surface  is  generated  as  a  copy  of  the  surface  of  the  gear  head-cutter.  Usually,  straight  profile 
blades  are  applied  for  the  gear  head-cutter.  Such  blades  form  the  head-cutter  generating  surface  as  a  cone. 

Application  of  parabolic  profiles  of  the  blades  instead  of  straight  line  profiles  permits  modification  of  the 
gear  tooth  surface  in  favor  of  conjugation  of  pinion-gear  tooth  surfaces. 

Simulation  of  Meshing.  The  algorithm  of  computerized  simulation  of  meshing  [8,  9]  provides  continuous 
tangency  of  contacting  gear  tooth  surfaces  Xj  and  E2  as  equality  of  position  vectors  and  unit 

normals  =  nf  ^ .  Instead  of  equality  of  unit  normals  to  contacting  surfaces,  collinearity  of  surface 

normals  may  be  considered.  Position  vectors  (i=l,2)  and  unit  normals  {i-h2)  are  considered  in 
fixed  coordinate  system  (Fig.  1). 

Application  of  TCA  determines  the  path  of  contact  on  the  pinion-gear  tooth  surfaces,  dimensions  and 
orientation  of  instantaneous  contact  ellipse  at  each  point  of  contact,  the  transmission  function 
relates  the  angles  of  rotation  of  the  gear  and  the  pinion,  and  function  A02(0i)  of  transmission  errors  is 
determined  as 


where  and  N2  are  the  numbers  of  teeth  of  the  pinion  and  the  gear. 

Predesign  of  Parabolic  Function  of  Transmission  Errors  [8].  The  output  of  TCA  for  a  misaligned  gear 
drive  shows  that  the  transmission  function  ^2(^1)  ^  piecewise  discontinuous  almost  linear  function  (Fig. 

2(a)).  The  transfer  of  meshing  with  the  transmission  function  ^2(0i)  (Pig-  2(^))  function  A^2(^i)  of 
transmission  errors  (Fig.  2(b))  is  accompanied  with  the  jump  of  the  angular  velocity  and  therefore  noise  and 
vibration  are  inevitable. 

The  purpose  of  synthesis  of  a  low-noise  gear  drive  is  to  transform  the  transmission  function  02  (0i) 
shown  in  Fig.  3(a)  and  obtain  a  parabolic  function  of  transmission  errors  shown  in  Fig.  3(b).  The  modified 
function  A02(0i)  (Fig.  3(b))  is  a  continuous  one  and  is  determined  as 

2 

(^i)  =  -«0i  (2) 


where  a  is  the  parabola  coefficient. 

It  is  important  to  recognize  that  the  function  of  transmission  errors  must  be  negative  and  therefore  the  gear 
will  lag  with  respect  to  the  pinion.  Therefore,  the  transfer  of  meshing  is  accompanied  with  elastic 
deformations,  and  as  a  result,  the  gear  ratio  is  increased. 

The  desired  parabolic  function  of  transmission  errors  is  provided  by  the  respective  synthesis  of  the  gear 
drive.  The  idea  of  such  synthesis  is  based  on  the  predesign  of  a  parabolic  function  obtained  by  application 
of  modified  roll  during  the  process  of  generation.  Modified  roll  is  provided  by  application  of  nonlinear 
relations  between  the  angles  of  rotation  of  the  cradle  of  the  generation  machine  and  the  pinion  during  the 
process  of  pinion  generation. 

The  advantage  of  the  predesign  of  a  parabolic  function  of  transmission  errors  is  that  such  a  function  is  able 
to  absorb  the  linear  discontinuous  function  of  transmission  errors  caused  by  misalignment.  This  is 
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(2) 

illustrated  by  drawings  of  Fig.  4(a)  that  shows  a  linear  function  A^,  4(a)).  Function 

A(j>2\(pi)  is  predesigned  and  function  A<j)f\<t>f)  is  caused  by  misalignment.  The  sum  of  functions 
and  A^f  (0,)  is  represented  in  Fig.  4(a)  as  A\i/2(w,)  =  ■  The  parabola  coefficient  a  is  the 

same  in  functions  • 

(2) 

The  figure  confirms  that:  (i)  the  predesigned  parabolic  function  absorbs  indeed  linear  function  A02  (0i) » 

and  (ii)  the  resulting  function  slightly  dislocated  with  respect  to  origin  of  A02^\0i) ,  function 

Ay/^iy/^)  is  an  unsymmetrical  one.  Parameters  c  and  d  of  the  origin  of  Ava2(V^i)  determined  as 

d  -b~ lAa  and  c-blla  .  In  the  process  of  meshing,  whereas  several  cycles  of  meshing  are  considered 
(Fig.  4(b)),  the  resulting  function  Axf/2{y/i)  of  transmission  errors  becomes  a  symmetric  parabolic  function. 

The  authors  have  developed  computerized  synthesis  of  formate-cut  spiral  bevel  gears  that  is  based  on 
simultaneous  application  of  local  synthesis  and  simulation  of  meshing  and  contact  of  gear  tooth  surfaces. 

The  purpose  of  local  synthesis  is  to  determine  the  pinion  machine-tool  settings  considering  as  given:  (i)  the 
gear  machine-tool  settings,  and  (ii)  the  conditions  of  meshing  and  contact  at  the  mean  contact  point  M  such 
as  the  tangent  to  the  path  of  contact  at  M,  the  major  axis  of  the  contact  ellipse  at  M,  and  the  derivative  of 
the  function  of  transmission  errors. 

The  purpose  of  simulation  of  meshing  and  contact  is  the  determination  of  the  bearing  contact  and  function 
of  transmission  errors  knowing  the  pinion  and  gear  machine-tool  settings. 

The  applied  computational  procedures  are  iterative  processes  with  the  goals  to  obtain:  (i)  a  longitudinally 
directed  bearing  contact,  (ii)  a  parabolic  function  of  transmission  errors  with  limited  value  of  maximal 
errors,  and  (iii)  reduced  sensitivity  of  the  gear  drive  to  errors  of  alignment. 


3.  Derivation  of  Gear  Tooth  Surface  [7,  9] 

Applied  Coordinate  Systems:  Coordinate  system  is  the  fixed  one  and  it  is  rigidly  connected  to  the 

cutting  machine  (Fig.  5),  Fig.  5(a)  and  5(b)  show  two  sets  of  coordinate  systems  applied  for  generations  of 
left  and  right  hand  gears,  respectively.  Coordinate  system  ^2  is  rigidly  connected  to  the  gear.  Coordinate 

system  is  rigidly  connected  to  the  gear  head-cutter.  It  is  considered  that  the  gear  head-cutter  is  a  surface 
of  revolution,  and  the  rotation  of  the  head  cutter  about  the  X^  -  axis  does  not  affect  the  process  of 
generation. 

The  installment  of  the  tool  on  the  machine  is  determined  by  parameters  H2  and  V2 »  called 

horizontal  and  vertical  settings.  Parameters  and  represent  the  settings  of  the  gear. 

Head  Cutter  Surface:  The  blades  of  the  gear  head-cutter  of  a  parabolic  profiles  are  shown  in  Fig.  6(a). 
Each  side  of  the  blade  generates  two  sub-surfaces.  The  segment  a  of  the  parabolic  profile  generates  the 
working  part  of  the  gear  tooth  surface.  The  circular  arc  of  radius  generates  the  fillet  of  the  gear  tooth 
surface.  The  generating  surfaces  of  the  head-cutter  are  formed  by  rotation  of  the  blade  about  the  X^  -  axis 
of  the  head-cutter  (Figs.  6(b)  and  6(c));  the  rotation  angle  is  9^  .  Therefore,  the  generating  surfaces  are:  (i) 

surface  of  revolution  formed  by  rotation  of  the  blade  of  parabolic  profile,  and  (ii)  the  surface  of  the  torus 
formed  by  rotation  of  the  circular  arc  profiles.  A  point  on  the  generating  surface  is  determined  by 
parameters  and  9^  for  the  working  surface,  and  by  and  9^  for  the  fillet  surface.  Angle  is  formed 
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between  the  tangent  line  of  the  blade  at  point  M  and  the  vertical  center  line  of  the  blade.  Parameter 
measured  from  point  M  in  the  chosen  direction  is  considered  as  a  positive  one,  and  angles  and  as 

the  acute  ones.  The  apex  of  the  parabola  is  located  at  point  M  determined  by  parameter  . 

The  surface  of  revolution  and  the  torus  surface  of  the  head-cutter  are  designated  as  parts  (a)  and  (b)  oi  the 
head-cutter  generating  surface.  Surface  of  the  head-cutter  is  represented  by  vector  function 

T^\ug,9g)  as  follows 


2  . 

-u,cosa^-\-aM^sma-s^  cosa^ 

6  5^6  <5  60  6 

{Rg  ±  sin  ±  sin  ±  cos  )  sin  9^ 

{Rg±  sin  ±  sin  ±  cos  )  cos  9^ 


(3) 


where  and  0,  are  the  surface  coordinates,  a,  is  the  blade  angle  at  point  M,  R^  is  the  cutter  point 

So  o  6 

radius  (Fig.  6)  and  given  by 

(4) 

The  upper  and  lower  signs  in  equations  (3)  and  (4)  correspond  to  the  concave  and  convex  sides  of  the  gear 
tooth,  respectively. 

The  unit  normal  to  the  gear  generating  surface  is  represented  by  the  equations 


nfiu^,9J- 


N 


(.0) 


N 


(a) 


(a) 


(a) 


X- 


du„  d0„ 


s 


(5) 


Equations  (3)  and  (5)  yield 


^^l  +  4ay^  (6) 


Surface  is  represented  in  as 


±  (sina^  +  la^Ug  cosa^ ) 
(cos  -  la^Ug  sin  )  sin  6^ 

(cos  -  2a^Ug  sin  )  cos  0^ 


rTiW 


-p„(l-cos\) 

(X^±p„sinAjsin0^ 

(Z^±p„sinA„)cos0^ 


(7) 


where 


Xg=Rg+  p^(l  -  sina^  )/cosa^ 
Here,  p^  is  the  fillet  radius  of  the  gear  generating  surface. 
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The  unit  normal  to  the  gear  generating  surface  is  represented  by  the  equations 


n  (X  e  )  = 

IN’ 


ib) 


(b) 


(b) 


N 


■(b) 


dr. 


(b) 


dK,  d0„ 


(8) 


Equations  (8)  yield 


nTaM= 


±  cosA^,, 
sinA^^,  sin^^ 
sinA„.  cos0^ 


(9) 


Equations  of  Head-Cutter  Surfaces  in  52:  Applying  coordinate  transformation  from  5^  to  52,  the  gear 
tooth  surfaces  can  be  represented  in  52  by  the  equations 


Here 


M,,  =M,„,  M„,^^ 


M 


10  0  0 

0  1  0  ±y, 

0  0  1  //j 

0  0  0  1 


(10) 

(11) 


(12) 


’cosy„,^  0  -siny„,^  0 

0  10  0 

M2„,  = 

^  sin  7,^  0  cosy^^ 

0  0  0  1 

Here,:  V2 ,  H2  and  are  the  gear  machine  tool  settings.  The  upper  and  lower  signs  in  front  of  V2 
correspond  to  right  hand  and  left  hand  gears  respectively.  The  whole  set  of  gear  machine  tool  settings  is 
represented  in  Table  1 .  Remember  that  the  generated  gear  tooth  surface  is  a  copy  of  the  surface  of  the  head- 
cutter,  that  is  a  surface  of  revolution.  The  rotation  of  the  head-cutter  about  the  X^  —  axis  that  is  necessary 
for  the  cutting  or  grinding  process  does  not  affect  the  shape  of  gear  tooth  surfaces. 

The  geometric  model  of  gear  tooth  surfaces  is  shown  in  Fig.  7. 

Table  1:  Gear  Machine  Tool  Settings 


Blade  angle  (Fig.  6) 
Blade  parabola  coefficient 
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Cutter  radius  (mm)  (Fig.  6) 

Point  width  (mm)  (Fig.  6) 

Cutter  point  radius  (mm)  (Fig.  6) 
Horizontal  settings  //2  mm)  (Fig.  5) 
Vertical  settings  V2  (Fig*  5) 

Machine  center  to  back  Xq  (mm)  (Fig.  5) 
Machine  root  angle  Jm,  (Fig-  5) 

Fillet  radius  (mm)  (Fig.  6) 


4.  Derivation  of  Pinion  Tooth  Surface  [7,  9, 10] 

Applied  Coordinate  Systems:  Coordinate  systems  applied  for  generation  of  pinion  are  shown  in  Fig.  8 
and  9.  Fig.  9(a)  and  (b)  correspond  to  generation  of  right  and  left  hand  pinions,  respectively.  Coordinate 
systems  and  are  the  fixed  ones  and  they  are  rigidly  connected  to  the  cutting  machine.  The 

movable  coordinate  systems  5j  and  are  rigidly  connected  to  the  pinion  and  the  cradle,  respectively. 
They  are  rotated  about  the  -  axis  and  -  axis,  respectively,  and  their  rotations  are  related  with  a 
polynomial  function  if  modified  roll  is  applied  (see  below).  The  ratio  of  instantaneous  angular 

velocities  of  the  pinion  and  the  cradle  is  defined  as  ))  =  The  magnitude  of 

¥ci  =0  is  called  ratio  of  roll.  Parameters  ,  X^^ ,  ,  and  are  the  basic  machine  tool 

settings  for  pinion  generation. 

Head  Cutter  Surfaces:  The  pinion  generating  surfaces  are  formed  by  surface  and  generated  by 

straight-line  and  circular  arc  parts  of  the  blades  (Fig.  10).  Surface  is  represented  as 


r:\sM 


(/?p  +  SpSma^)cos0p 
(/?pTSpSino:p)sin0p 
-SpCosa^ 


(13) 


where  and  6^  are  the  surface  coordinates;  is  the  blade  angle;  is  the  cutter  point  radius  (Fig.  10). 

The  upper  and  lower  signs  in  equations  (13)  correspond  to  the  convex  and  concave  sides  of  the  pinion  tooth 
respectively. 

The  unit  normal  to  the  pinion  generating  surface  '  is  represented  by  the  equations 


nf{ep)  = 


N' 


(a) 


N' 


(a) 


dr| 


(a) 


Equations  (13)  and  (14)  yield 


n;'(ej  = 


cosa„  cos^„ 

p  p 

cosa^  sin^^ 
Tsina^ 


(14) 


(15) 


NASA/CR— 2001-210894 


8 


For  the  fillet  surface  we  obtain 


(X^  +  pj  sinAy.)cos^p 


{Xj,e^)=  (X^ +  p^sinA^)sin0p 
-  Pj{\-cosXj) 


where 


X  ^-Rp±  p^(l  +  sina^)/cosa^ 


and  p.  is  the  radius  of  the  tool  fillet. 


The  unit  normal  to  the  pinion  generating  surface  2  is  represented  by  the  equations 


(b^  (b) 


Equations  (17)  yield 


Fsin  cos0^ 


Hp  (^j,0p)=  sin A^ sin (18) 
+  cosA^ 

The  head-cutter  is  mounted  on  coordinate  system  5^^  called  the  cradle  of  the  cutting  machine  and  its 
installment  is  determined  by  settings  5,.  and  (Fig.  9).  In  the  process  of  pinion  generation  coordinate 
systems  and  performed  related  rotations  with  respect  to  .  Then,  in  coordinate  system  5j ,  a 
family  of  pinion  head-cutter  surfaces  is  generated  and  determined  as: 

(19) 

(20) 

where  y/^  is  the  generalized  parameter  of  motion,  and 


10  0  5,  cosg, 

0  10  ±S.  sin^i 

M  =  ' 

0  0  1  0 

0  0  0  1 
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COSI/^^^ 

-sinyr^, 

0 

0 

sinv^,, 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

M 


1  0 
0  1 
0  0 
0  0 


0 

0 

1 

0 


0 


1 


M 


M,,. 


siny„,, 

0 

-cosy„,, 

0 

0 

1 

0 

0 

cos7„„ 

0 

siny,,,, 

- 

0 

0 

0 

1 

r  cosv^. 

sini/r,  0 

“I 

-sinv^i  cosv/'i  0  0 
0  0  10 
0  0  0  1 


The  upper  and  lower  signs  of  5^  sin^i  in  matrix  correspond  to  generation  of  right  and  left  hand 
pinions,  respectively.  Where  the  modified  roll  is  applied  in  the  process  of  generation,  the  rotation  angles 
i/Aj  and  of  the  pinion  and  cradle  are  related  as 


Wi  +b2wl  +hwl 


(21) 


where  ^2  modified  roll  coefficients. 

The  derivative  of  function  )  taken  at  =0  determines  the  so-called  ratio  of  roll  represented  in 

equation  (21)  by  . 


Equation  of  Meshing:  The  pinion  tooth  surface  Zj  is  the  envelope  to  the  family  of  head-cutter  surfaces. 
Surface  Xj  consists  of  two  surfaces  and  x|^^  which  correspond  to  the  cutter  surfaces  and  , 

respectively.  Modified  roll  is  applied  in  the  process  of  generation.  For  generation  of  surface  Xj"*^ ,  the 
equation  of  meshing  is  represented  as 


(a)  (pi)  \  n 

“m, =flp(Sp,e^,y/^)  =  0 


(22) 


where  the  unit  normal  to  the  surface,  and  is  the  velocity  in  relative  motion.  These  vectors  are 
represented  in  the  fixed  coordinate  system  as  follows: 


(a)  y 

n,  =L 

nil 


miC,  ^CiP‘ 


n:\ej 


(23) 
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(24) 


Here 


(1)^ 


vir  =[(«ir-w™’)xr„,,  -(0,„  xw';;) 


r,, 

o„,o,,  =[o.  X«J  ' 

®I!*  =  [cosr,„,.  0,  sin7„J 


The  ratio  w^j  is  not  constant  since  modified  roll  is  applied  and  can  be  represented  as 

ft),,  d\(f  /dy/,  1 

=— !-  =  — -i-/— —  = - j 

(Oi  at  /  dt  mi^  +  2Z72V^i-f3Z?3V^i 

where  b2  and  are  modified  roll  coefficients. 

Finally,  we  obtained  the  equations  for  pinion  tooth  surface  part  (a)  as 

Similarly,  the  fillet  surface  of  pinion  tooth  surface  may  be  represented  as 

,6^) 

/;;\A,,0^,VA,)  =  O 

The  geometric  model  of  pinion  tooth  surfaces  is  shown  in  Fig.  11. 


(25) 


(26) 

(27) 


(28) 

(29) 


5.  Procedures  of  Proposed  Design 

The  authors  have  developed  an  approach  for  the  design  of  formate  cut  spiral  bevel  gears  that  provides  a 
stabilized  bearing  contact  and  reduced  level  of  noise.  The  bearing  contact  is  of  a  longitudinal  direction  that 
enables  to  increase  the  contact  ratio,  avoid  edge  contact  and  reduce  the  shift  of  the  bearing  contact  caused 
by  misalignment.  The  reduction  of  noise  is  achieved  by  application  of  a  predesigned  parabolic  function  of 
transmission  errors  and  limitation  of  maximal  transmission  errors  caused  by  misalignment.  The  approach  is 
based  on  three  procedures  that  require  simultaneous  execution  of  computer  programs  developed  for  the 
local  synthesis  and  simulation  of  meshing  and  contact  of  the  gear  drive. 

The  purpose  of  the  local  synthesis  is  to  determine  the  pinion  machine-tool  settings  for  the  following 
conditions  of  meshing  and  contact  of  pinion-gear  tooth  surfaces  Zj  and  Z2 : 
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(i)  The  gear  machine-tool  settings  are  considered  as  given  (they  may  be  adapted,  for  instance,  from  a 
Gleason  summary  for  gear  generation); 

(ii)  The  pinion-gear  tooth  surfaces  Zj  and  X2  are  in  tangency  at  a  chosen  mean  contact  point  M  (Fig.  12); 

(iii)  Parameters  2a,  %  ,  and  m[2  are  assigned  ahead  and  determine  the  conditions  of  meshing  and  contact  of 
E,  and  X2  at  point  M  and  in  the  neighborhood  of  M  Here,  2a  is  the  major  axis  of  the  instantaneous  contact 

ellipse,  7J2  determines  the  tangent  to  the  path  of  contact  on  the  gear  tooth  surface  at  M,  and  m[2  is  the 
derivative - at  point  M,  where  rj\2((j)i)  is  the  gear  ratio. 

Note:  It  will  be  shown  in  section  6  that  the  location  of  point  M2  on  the  gear  tooth  surface  as  the  candidate 
for  the  mean  contact  point  of  Xj  and  Z2  may  be  chosen.  Then,  the  procedure  of  the  local  synthesis  permits 
the  determination  of  point  Mj  on  pinion  tooth  surface  that  will  be  in  tangency  with  Z2  the  chosen 
point  M2  of  X2  • 

The  simulation  of  meshing  and  contact  of  pinion-gear  tooth  surfaces  is  accomplished  by  application  of  the 
TCA  (Tooth  Contact  Analysis)  computer  program  developed.  The  input  for  TCA  are  the  equations  of 
pinion-gear  tooth  surfaces,  parameters  of  motion  and  assembly.  The  output  is  the  bearing  contact  and  the 
function  of  transmission  errors. 

The  design  is  based  on  simultaneous  execution  of  computer  programs  developed  for  local  synthesis  and 
TCA,  it  is  an  iterative  process  and  requires  application  of  following  three  procedures. 

Procedure  1: 

The  path  of  contact  on  gear  tooth  surface  is  a  spatial  curve  L.  Projection  of  L  on  tangent  plane  T  (Fig.  13)  is 
designated  by  Lj. .  The  purpose  of  procedure  1  is  to  obtain  that  is  a  straight  line  directed  longitudinally. 

Fig.  13  shows  projections  and  that  deviate  from  the  desired  line  Lj- . 

Fig.  14  shows  the  flow  chart  that  describes  the  steps  of  procedure  1.  The  procedure  of  computation  is  an 
iterative  process  and  requires  variation  of  parameter  m\2 .  The  steps  of  the  procedure  are  as  follows: 

Step  1:  Input  the  gear  machine-tool  settings  into  the  program  of  local  synthesis. 

Step  2:  Input  parameters  2a,  1)2 ,  and  m^2  ^^^o  the  program  of  local  synthesis. 

Step  3:  Execute  the  program  of  local  synthesis. 

Step  4:  Determine  from  the  program  of  local  synthesis  the  pinion  machine-tool  settings. 

Step  5:  Using  the  machine-tool  settings  obtained,  compute  the  pinion  and  gear  tooth  surfaces  Xj  and  X2 . 
Step  6:  Execute  TCA  computer  program  and  obtain  projection  of  path  of  contact  L  on  tangent  plane  T 
(Fig.  13). 

Step  7:  Execute  the  iterative  process  and  obtain  projection  of  path  of  contact  L^as  a  straight  line  by 
variation  of  m^2  • 

Step  8:  Output  the  final  machine-tool  settings  for  the  pinion. 

The  iterative  process  is  based  on  the  following  considerations: 

(1)  The  TCA  computer  program  applied  determines  projection  of  path  of  contact  numerically. 

(2)  Consider  coordinate  system  (Fig.  13)  which  origin  coincides  with  M  and  axes  (x^y,)  are 
located  in  plane  T .  Axis  x,  is  the  tangent  to  at  point  M  . 
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(3)  Represent  projection  as  a  parabolic  regression  equation  by  using  a  subroutine  of  regression  [11] 


(30) 


where  m\^  is  the  input  parameter  that  is  varied  in  the  process  of  iteration 

(4)  The  goal  is  to  represent  as  a  straight  line  that  might  be  obtained  by  such  a  magnitude  of  m\2 

whereas  coefficient  (see  equation  (30)).  The  iterations  for  determination  of  =  0  are 

executed  by  applying  the  secant  method  [12]  and  illustrated  by  Fig.  15.  Equation  (30)  with  ^2(^12)“^ 
provides  the  projection  of  path  of  contact  as  a  straight  line. 


Procedure  2: 

Procedure  1  discussed  above  provides  the  desired  longitudinal  path  of  contact  and  bearing  contact. 
However,  the  shape  of  the  obtained  function  of  transmission  errors  and  the  magnitude  of  maximal 
transmission  errors  do  not  satisfy  the  requirements  of  low-noise  gear  drive.  The  goal  of  procedure  2  is  to 
obtain  a  parabolic  function  of  negative  transmission  errors  and  of  limited  value  of  maximal  errors. 

The  sequence  of  steps  applied  for  procedure  2  is  represented  by  the  flow  chart  shown  in  Fig.  16. 

At  the  start,  we  use  the  pinion  and  gear  machine  tool  settings  obtained  in  procedure  1,  but  consider  that  the 
pinion  is  generated  by  application  of  modified  roll.  Then,  we  compute  the  gear  and  pinion  tooth  surfaces 
Zj  and  Z2 ,  and  use  them  for  the  TCA  computer  program.  The  intermediate  output  of  TCA  is  a  function  of 
transmission  errors.  We  compare  it  with  the  sought-for  parabolic  function  of  transmission  errors  and 
determine  the  necessary  corrections  using  modified  roll.  Repeating  the  procedure  discussed  above,  we  can 
obtain  the  sought-for  parabolic  function  of  transmission  errors. 

Analytically,  the  algorithm  of  procedure  2  is  represented  as  follows. 

Step  1:  Function  of  transmission  errors  A02(0i)  obtained  from  TCA  is  numerically  represented  as  a 
polynomial  function  up  to  third  order  to  be  included  [11] 


(31) 


The  shape  and  magnitude  of  maximal  errors  of  function  A02  (A)  satisfy  the  requirements  of  a  low- 

noise  gear  drive.  Function  (31)  has  to  be  transformed  into  a  parabolic  function  with  limited  magmtude  of 
transmission  errors  by  application  of  procedure  2. 

Step  2:  We  apply  for  transformation  of  function  A^^^^j)  the  modified  roll  for  generation  of  the  pinion 
that  is  represented  as 


Wi  (¥c, )  =  nh,¥c,  +  hwl  +  hwl  (32) 

Here,  i/r,  and  \ff2  are  the  angles  of  rotation  of  the  pinion  and  the  cradle  performed  during  the  process  for 
generation,  is  the  first  derivative  of  taken  at  =0  and  obtained  by  the  procedure  of  local 
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synthesis  (see  section  6).  Coefficients  and  have  to  be  determined  by  using  an  iterative  process  (see 
below). 

Step  3:  Assigning  coefficients  Z?2  ^3  modified  roll  and  using  pinion  machine-tool  settings  obtained 

in  procedure  1,  we  may  determine  the  pinion  tooth  surface  Ej .  The  gear  tooth  surface  E2  is  considered  as 
given. 

Step  4:  Applying  TCA  for  simulation  of  meshing  of  Xj  and  X2  we  may  obtain  numerically  the  modified 
function  of  transmission  errors  and  then  represent  it  again  by  polynomial  function  (31)  with  new 
magnitudes  of  coefficients  Uq,  ^2  and  . 

Procedure  2  with  steps  1  to  4  has  to  be  repeated  until  a  parabolic  function  of  transmission  errors  A<^2(^i) 
with  limited  magnitude  of  maximal  transmission  errors  is  obtained. 

The  computation  is  an  iterative  process  based  on  application  of  secant  method  of  numerical  computations 
[12].  The  goal  is  to  transform  the  shape  of  obtained  function  A<l>2\(l>i)  of  transmission  errors  into  a 
parabolic  function 


A02(0i)  =  -a20i  , 


■— <0j<  — 


where 


|a02(^,)| 


max 


=  AO 


(33) 


(34) 


This  goal  is  accomplished  by  variation  of  coefficients  Z?2  and  Z73  of  function  (32)  provided  by  modified 
roll.  The  variation  of  Z?2  h  performed  independently  and  is  illustrated  by  drawings  of  Fig.  17. 


Fig.  17(a)  illustrates  variation  of  coefficient  b^  of  modified  roll  for  obtaining  of  coefficient  ^3—0  of 
function  (31).  Function  a2(^2)  determined  as  the  output  of  TCA  by  variation  of  modified  roll. 

Fig.  17(b)  illustrates  variation  of  coefficient  Z?2  of  modified  roll  for  obtaining  of  |A02(0i)|roax  “  * 


Function  A02(0i)  is  determined  as  the  output  of  TCA  by  variation  of  Z?2 . 


Procedure  3: 

Procedures  1  and  2  are  performed  by  simultaneous  application  of  computer  programs  of  local  synthesis  and 
TCA,  but  the  errors  of  alignment  of  the  gear  drive  are  not  considered. 

The  purpose  of  procedure  3  is  to  adjust  the  gear  drive  for  the  existence  of  errors  of  alignment.  This  is 
achieved  by  adjusting  the  previously  designed  path  of  contact  to  the  assigned  or  expected  errors  of 
alignment.  The  procedure  is  based  on  the  flow  chart  shown  in  Fig.  18. 

The  adjustment  of  the  gear  drive  to  the  errors  of  alignment  is  performed  by  correction  of  parameter  772  (Fig. 
12)  that  determines  the  orientation  of  the  path  of  contact  on  the  gear  tooth  surface. 
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The  prescribed  procedure  yields  that  in  some  cases  of  design  it  becomes  necessary  to  deviate  the  path  of 
contact  from  the  longitudinal  direction  to  reduce  the  shift  of  bearing  contact  by  errors  of  alignment.  The 
procedure  has  to  be  applied  for  separate  simulation  of  each  error  of  alignment  for  obtaining  the  required 
correction  of  rj,  •  The  permissible  tolerance  of  each  alignment  error  is  obtained  by  procedure  3. 


Local  Synthesis  and  Determination  of  Pinion  Machine  Tool  Settings  [7,  8, 9] 

Remembering  that  the  local  synthesis  is  accomplished  simultaneously  with  the  TCA,  and  the  computation 
is  an  iterative  process.  The  computer  program  developed  permits  the  determination  at  each  iteration  the 
pinion  machine-tool  settings. 

The  input  parameters  for  the  local  synthesis  are  the  gear  machine-tool  settings  and  parameters  rj^  ,  a  and 
mji  (see  Flow  chart,  Fig.  14). 

The  procedure  of  the  local  synthesis  is  as  follows. 

Stage  1:  Determination  of  the  mean  contact  point  on  the  gear  tooth  surface 

Step  1:  Mean  point  A  on  surface  Ej  is  chosen  by  designation  of  parameters  and  (Fig.  19),  where 
A  is  the  candidate  for  the  mean  contact  point  M  of  surfaces  Ej  and  Ej .  Then  we  obtain  the  following 
system  of  two  equations  in  two  unknowns 


X2(M*,F*)  +  y,'(u’,F*)  =  R' 


(35) 


where  Xj,  Tj  ^nd  are  the  projections  of  position  vector 

*  * 

Step  2:  Equations  (35)  considered  simultaneously  permits  the  surface  parameters  {Ug,9^)  to  be  determined 

for  point  A.  Vector  functions  r^(Ug,6^)  and  n^{dg)  determine  the  position  vector  and  surface  unit  normal 

*  * 

for  a  current  point  of  gear  tooth  surface  .  Taking  in  these  vector  functions  and  6^  —9^,  we  can 

determine  the  position  vector  of  point  A  and  the  surface  unit  normal  at  A. 

Step  3:  Parameters  m*  and  9^  and  the  unit  vectors  and  of  principal  directions  on  surface  are 
considered  as  known.  Here: 


e 


8 


-V  (tj) 

-  cosa^  +  lau.  sina„ 

dUg 

- 

(±  sin  ag  ±  2a ^Ug  cosa^ )  sin  9^ 

(±sina^  ±  2a^Ug  cosa^  )cos0j 

(36) 
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e, 


(37) 


(a) 

■ 


-N  {(f) 

or. 

0 

s 

= 

cos^^ 

do. 

8 

s 

-sm0^_ 

The  gear  generating  surface  is  a  surface  of  revolution,  and  the  principal  curvatures  and  of  are 
determined  by  the  following  equations: 


= 


2a. 


3 

(l  +  4alulf 


-  cosa^  +  2a  sina^ 


(38) 


{Rg  ±  sin  ±  sin  ±  a^u^  cos  )yl-h4^^ 


The  upper  and  lower  signs  correspond  to  the  concave  and  convex  sides  of  gear  tooth,  respectively. 

The  unit  vectors  and  of  principal  directions  on  surface  S2  obtained  by  using  matrix 
transformation  from  to  82-  The  gear  principal  curvatures  and  of  E2  are  the  same  as  the 
principal  curvatures  and  ,  respectively,  of  the  gear  head-cutter. 

Stage  2:  Tangency  of  Surfaces  X2  (2^^  )  and  Xj  at  Mean  Contact  Point  M, 

(/i) 

Step  1:  The  derivations  accomplished  at  stage  1  permits  the  position  vector  T2  to  be  determined  and  the 

surface  unit  normal  112^^  of  point  A  of  tangency  of  surfaces  X2  (X^  ).  The  goal  now  is  to  determine  such  a 
point  M  in  the  fixed  coordinate  system  5^  (Fig.  20)  where  two  surfaces  X2  (X^ )  and  X^  are  in  tangency 
with  each  other. 

Surfaces  X2  and  X^  are  the  same  ones,  because  there  is  no  relative  motion  between  the  cradle  and 
workpiece  (between  coordinate  system  82  to  8^  whereas  the  gear  is  generated.  Using  coordinate 

transformation  from  52  to  5/  (Fig.  20),  we  may  obtain  and  .  The  new  position  of  point  A  in  5^ 
will  become  point  M  of  tangency  of  X2  and  Xj ,  if  the  following  equation  of  meshing  between  X2  and  Xj 
is  observed 


(A) 


,(0) 


(2UA)^(0) 


(02  )*V/  ’  (0 


)  =  0 


(39) 


Here,  and  ^  ^  ^(12, a)  relative  velocity  at  point  A  determined  with  the  ideal 

gear  ratio 


(0) 

AM21  = 


(2) 


0) 


(1) 


0) 


(40) 


The  solution  of  equation  (39)  for  02^^  provides  the  value  of  the  initial  turning  angle  for  02  •  It  is  evident 
that  two  surfaces  X2  and  X^  are  now  in  tangency  with  each  other  at  point  M. 
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Step  2:  We  consider  as  known  at  point  M  the  principal  curvatures  and  of  surface  E2 ,  and  the  unit 
vectors  and  of  principal  directions  on  Zo .  The  unit  vectors  and  are  represented  in  5^ .  The 
goal  now  is  to  determine  at  point  M  the  principal  curvatures  and  kj  and  kj^  of  surface  ,  and  the  unit 
vectors  and  of  principal  directions  on  Xj .  This  goal  can  be  achieved  by  application  of  the 
procedure  described  in  Appendix  1 .  It  is  shown  in  Appendix  1  that  the  determination  of  e^-  and 

C/,  becomes  possible  if  parameters  ,  r}^  (or  77, ),  and  a! 5  are  assumed  to  be  known  or  are  used  as 
input  data. 

Stage  3:  Tangency  of  Surfaces  X2  ( ) ,  Xj  and  X^  at  Mean  Contact  Point  M. 


Basic  Equations:  Tangency  of  X2  (X^ )  and  Xj  at  mean  contact  point  M  has  been  already  provided  at 

the  previous  stages.  The  position  vector  of  point  M  and  the  surface  unit  normal  at  point  M  were 
determined  in  coordinate  system  .  Let  us  imagine  now  that  coordinate  system  Sj  that  coincides  with 

and  surface  Xj  are  installed  in  coordinate  system  (Fig.  8).  Let  angle  be  the  installment  angle  of 
the  pinion  (it  is  the  initial  value  of  y/^ ).  Using  coordinate  transformation  from  to  ,  we  may 

determine  in  position  vector  of  point  M  and  the  surface  unit  normal  .  In  Appendix  1  the 
conditions  of  improved  meshing  and  contact  of  pinion  and  gear  tooth  surfaces  X^  and  X2  are  considered 
and  then  the  relationships  between  the  principal  curvatures  and  directions  of  surfaces  for  such  conditions  of 
meshing  and  contact  are  determined  (see  equations  (93)).  The  point  of  tangency  of  surfaces  Xj  and  X^  is 

designated  in  Appendix  1  as  point  B.  The  pinion  generating  surface  X^  is  installed  in  taking  the 
cradle  angle  y/^  equal  to  zero.  The  position  vector  of  point  B  of  surface  X^  and  the  surface  unit  normal  at 

B  are  represented  in  as  and  .  The  tangency  of  Xj  and  X^  at  the  mean  contact  point  M  is 
satisfied,  if  the  following  vector  equations  are  observed 


(A/)  (B) 

n.„  =  n„. 


(41) 


(M)  (B) 

r„.  =  r... 


(42) 


(M)  Up)  . 

n„.  •  =0 


(43) 


where  equation  (43)  is  the  equation  of  meshing.  Observation  of  equations  (41)-(43)  means  that  all  of  the 
three  surfaces  ( X2  ( X^ ) ,  X^  and  X^ )  are  in  contact  at  point  M. 

Using  equations  (41)-’(43)  and  (93)  from  Appendix  1,  it  becomes  possible  to  obtain  the  settings  of  the 
pinion  and  the  head-cutter  that  guarantee  the  improved  conditions  of  meshing  and  contact  at  point  M.  The 
machine  tool  settings  to  be  determined  are  as  follows: 


(i)  Xg,  and  Settings  and  X^^  are  related  by  the  following  equation  (Fig. 

21) 

=-(^z),-i^)sinr„„  (44) 


where 
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-[(A„  +^»V2)cosy,  sinyi] 


\o,o\  = 


[(A„  +  F.,/2)smyi-a^  cosy,] 
tany„,_ 


Equation  (44)  is  derived  for  the  case  of  design  of  a  spiral  bevel  gear  with  different  apexes  of  the  pitch  and 
root  cones. 

(ii)  Design  parameter  of  the  head-cutter  (Fig.  10). 

(iii)  Parameters  5,.  and  that  determine  the  installment  of  the  head-cutter  on  the  cradle  (Fig.  9). 


(iv)  Parameter  is  the  initial  angle  for  installment  of  coordinate  system  with  respect  to  5^^ 
(Fig.  8),  and  parameter  of  the  head-cutter  surface  . 

Computation  of  Pinion  Machine-Tool  Settings.  The  procedure  for  computation  is  as  follows: 

Step  1:  Determine  the  values  of  6^  and  that  are  the  sought-for  two  unknowns.  Equation  (41)  is  used 
for  determination  of  0^  and  ,  taking  into  account  that 


and 


cosa^  cosOp 
cosa^  sin0^ 
+  sin 


(M) 


j  y  (M) 

1^1 


where  nj^^  since  coincides  with  Sj  (Fig.  20).  Here, 


^m^b^ 


L,,= 


sin7„. 

0 

cosy,,, 

0 

1 

0 

-cosr^^ 

0 

siny„,. 

(0) 

COSlf/^ 

(0) 

(0)  . 

-siny/'i  ( 

_ (0)  , 

smi/A, 

0 


COS^/’l 
0 


(45) 


(46) 


(47) 


(48) 


Equations  (41)  and  (45)-(48)  yield  the  following  expressions  for  determination  of  6^  and 


(0) 


n^  +  siny  sina 
COS0„  = - - - - 

cosy  cosa 


(49) 


.  ,0)  -n,, cose  cosa  +  sin y  sin9p  cosa^  +  cosy„,^ 

suiyr,  —  2  2. 

(n,,  +  n,,,)siny„,^ 


(50) 
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(51) 


(0) 

cosy/ 1 


sin^  cos  a  -  sini/A, 


(0) 


where  is  the  given  value  of  the  profile  angle  of  the  head-cutter,  and  and  are  the  three 

components  of  vector  The  great  advantage  of  the  approach  developed  is  that  the  requirement  of  the 

coincidence  of  the  normals  does  not  require  a  non-standard  profile  angle  or  the  tilt  of  the  head-cutter 
with  respect  to  the  cradle. 

Using  Op  ,  it  becomes  possible  to  determine  the  unit  vectors  and  e,  of  principal  directions  on  surface 
Zp  at  point  M. 

Step  2:  Determination  of  pinion  machine-tool  settings  )»  ’  ^i^and  the  design  parameter  Rp 

of  the  head-cutter  (five  unknowns  are  sought-for). 

As  a  reminder,  and  are  related  by  equation  (44).  The  determination  of  the  machine-tool  settings 
mentioned  above  is  based  on  application  of  the  system  of  equations  (93)  and  equation  (43)  that  represent  a 
system  of  four  non-linear  equations  with  four  unknowns:  X^,^ ,  ,  m^p  and  Rp  . 

The  design  parameters  mentioned  above  provide  as  well  improved  conditions  of  meshing  and  contact  at  the 
mean  contact  point  M. 

Step  3:  Determination  of  machine-tool  settings  5^  and  q^  and  the  pinion  surface  parameter  Sp  (three 
unknowns  are  sought-for). 

Determination  of  the  mentioned  three  parameters  is  based  on  application  of  equation  (42),  considering  that 
generating  surface  is  a  cone.  The  final  equations  for  a  right  hand  pinion  are  as  follows: 


cosq,  +  (/?p  +  5p  sinap)cos0p  =  ’ 

(52) 

sin(±5,)  +  (i?p  +  5p  sinap)sinep  =  * 

(53) 

(M) 

—  S  COSOC  =  Z  1 

(54) 

The  upper  and  lower  signs  in  front  of  in  equation  (53)  correspond  to  the  design  of  right-hand  and  left 
hand  pinions,  respectively. 

The  stages  of  computation  discussed  above  may  be  summarized  as  follows: 

(i)  It  is  necessary  to  determine  ten  unknowns:  six  machine-tool  settings  ( X^  ,  ,  X^  ,  ,  m^p)  , 

two  surface  parameters  (Bp,  Sp ),  one  cutter  parameter  Rp  and  one  position  parameter  y/l^^  which  defines 
the  pinion  initial  turn  angle. 

(ii)  The  equation  system  for  determination  of  the  unknowns  is  formed  as  follows: 


(M)  (B) 

n„,  =  n„, 

nil 


(55) 
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II 

(56) 

(M)  (Ip)  r\ 

-V  =0 

(57) 

-(^D,  -  0^0,  )sin7„,_ 

(58) 

In  addition,  we  use  three  curvature  equations 


tan  2(7“"^= - 


k-k 


^23  ^13  ^/?)^33 


P  ,  . 

J33  sinztJ 


^13  ^"^3 


^33 


(59) 


Equation  (55)  is  equivalent  to  two  independent  scalar  equations,  equations  (56)  is  equivalent  to  three  scalar 
equations,  equations  (57),  (58)  and  (59)  represent  five  scalar  equations.  Thus,  the  system  of  equations 
provides  indeed  ten  scalar  equations  for  determination  of  ten  unknowns.  The  solution  for  the  unknowns 
requires  solution  of  a  subsystem  of  four  nonlinear  equations  (see  step  2)  and  solution  of  six  remaining 
equations  represented  in  echelon  form  (each  of  the  six  equations  contains  one  unknown  to  be  determined). 


7.  Simulation  of  Meshing  and  Contact  Accomplished  in  Procedure  1  (Section  5) 

[7,  8,  9] 

The  main  goal  of  simulation  of  meshing  and  contact  in  procedure  1  (see  Flow  chart  of  Fig.  14)  is  to 
determine  bearing  contact  that  corresponds  to  pinion  machine-tool  settings  obtained  in  section  6. 
Remembering  that  the  computational  procedures  represented  in  sections  6  and  7  must  be  applied 
simultaneously,  to  represent  an  iterative  process  for  determination  of  longitudinal  path  of  contact.  The 
developed  TCA  (Tooth  Contact  Analysis)  computer  program  (based  on  simulation  of  meshing  and  contact) 
provides  at  each  stage  of  iteration  the  path  of  contact.  In  addition,  TCA  provides  the  function  of 
transmission  errors  which  shape  has  to  be  corrected  by  application  of  pinion  modified  roll  (see  procedure  2 
of  section  5). 

Applied  Coordinate  Systems:  The  meshing  of  gear  tooth  surfaces  is  considered  in  the  fixed  coordinate 
system  that  is  rigidly  connected  to  the  housing  (Fig.  22).  Movable  coordinate  system  and  ^2  are 
rigidly  connected  to  the  pinion  and  the  gear,  respectively.  Auxiliary  coordinate  systems  and  5^^  are 
used  to  represent  the  rotation  of  the  pinion  (with  respect  to  5^^ )  and  the  gear  (with  respect  to  ).  The 
errors  of  alignment  are  simulated  by  respective  installment  of  5^^  and  with  respect  to  . 

The  errors  of  installment  are:  A//  —  the  axial  displacement  of  the  pinion;  Ay  —  the  change  of  the  shaft 
angle  7 ;  AE  -  the  shortest  distance  between  the  axes  of  the  pinion  and  the  gear  when  the  pinion-gear  axes 
are  crossed  instead  of  intersected;  AD  -  the  axial  displacement  of  the  gear.  In  the  case  of  an  aligned  gear 
drive  we  consider  that  A//  ,  Ay  ,  AD  and  AD  are  all  equal  to  zero. 
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Simulation  Algorithm:  During  the  process  of  meshing,  the  pinion  and  gear  tooth  surfaces  must  be  in 
continuous  tangency  and  this  can  be  provided,  if  at  any  instant  their  position  vectors  coincide  and  the 
surface  normals  are  collinear.  Instead  of  collinearity  of  surface  normals,  equality  of  surface  unit  normals 
may  be  required. 

Pinion  and  gear  tooth  surfaces  are  represented  in  coordinate  system  Sj,  by  the  following  equations 


,  (0,  )r, 


Ki  ’  ^2  )  ~  ^6, 2  (^2  )*'2  ) 


(60) 

(61) 


Here: 


M 


byi 


0 

0 


M,,,= 


COS</>2 
sin  02 
0 
0 


sin  01 

0 

0 

cos  01 

0 

0 

0 

1 

0 

0 

0 

1 

0  0 

0  ■ 

1  0 

0 

0  1  AH 

0  0 

Ij 

-sin  02 

0 

o' 

cos  02 

0 

0 

0 

1 

0 

0 

0 

1 

= 


cos(7  +  Ay)  0  -sin(7  +  A7)  -ADsin(7  +  A7) 
0  10  AE 

sin(7  +  A7)  0  cosCy  +  Ay)  ADcos(y  +  Ay) 
0  0  0  1 


Unit  normals  to  surfaces  of  pinion  and  gear  are  represented  in  S,^  by  the  following  equations,  respectively 


nr  (Sp  ’61  0, )  =  L,„  (^,  )n,  (5  6»  ) 


(2) 


=  L;,  2(02)n2(Mj,0^) 


(62) 

(63) 


Here: 


COS0J  sin0j  0 

-sin0j  COS0J  0 


0 


0  1 


NASA/CR— 200U210894 


21 


K2  = 


COS  02 

sin  02 

0 


-sin  02  0 
cos  02  0 

0  1 


L 


hbj 


cos(7  +  Ay)  0  -  sin(7  +  Ay) 
0  1  0 
sin(7  +  A7)  0  cos(7  +  A7) 


Conditions  of  continuous  tangency  of  pinion  and  gear  tooth  surfaces  are  represented  by  the  following 
equations 


(64) 

”^?(Sp,dp,Vi  A)-  A  ,<t>2)=^ 

(65) 

(66) 

Surfaces  Xj  and  X2  are  represented  in  Sf^  by  four  and  three  related  parameters,  respectively.  Equation 
(66)  is  the  equation  of  meshing  of  the  pinion  and  generating  head-cutter  and  relates  the  surface  parameters 
of  the  head-cutter  with  the  generalized  parameter  y/^  of  motion  of  the  process  of  generation  of  the  pinion. 
Vector  equations  (64)  ,  (65)  and  (66)  yield  three,  two  independent  and  one  scalar  equations,  respectively. 
The  whole  system  of  equations  (64)  -  (66)  provides  six  equations  for  determination  of  seven  unknowns 
represented  as 


/;.(j  ,V^i,0, 0^,02)  =  0,  fiGC,  (/  =  !,. ..,6) 


(67) 


One  of  the  unknowns,  say  0i ,  is  considered  as  the  input  parameter  in  the  range  - 


^  7C  ^ 


^<t>i 


N, 


.  It  is 


assumed  that  the  Jacobean  of  the  system  (67)  differs  from  zero  at  each  iteration. 

The  paths  of  contact  on  the  pinion  and  the  gear  tooth  surfaces  are  represented  by  the  following  functions: 

*1 


(68) 

(69) 


and 


^2  =r2(M^(^l),e^(0l)) 


(70) 


respectively. 

The  function  of  transmission  errors  is  defined  as 


N2 


(71) 
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The  bearing  contact  is  formed  as  a  set  of  instantaneous  contact  ellipses.  The  lengths  of  the  major  and  minor 
axes  of  contact  ellipse  and  their  orientation  are  determined  using  the  approach  proposed  in  [8]. 

Using  vector  function  (70),  we  may  determine  projection  of  the  path  of  contact  on  the  gear  tooth 
surface  on  the  tangent  plane  at  the  mean  contact  point  M.  The  goal  is  to  obtain  as  a  straight  line  directed 
longitudinally.  The  process  of  computation  is  an  iterative  process  and  requires  simultaneous  application  of 
local  synthesis  and  TCA.  The  errors  of  alignment  in  this  procedure  are  not  taken  into  account.  The 
influence  of  errors  of  alignment  on  transmission  errors  requires  additional  application  of  TCA, 

It  may  be  discovered  in  some  cases  of  design,  that  the  gear  drive  is  too  sensitive  to  errors  of  alignment. 
Then,  it  becomes  necessary  to  deviate  from  the  longitudinal  directed  straight  line  by  correction  rj,  in 
procedure  3  (see  Flow  chart  of  Fig.  18). 


8.  Numerical  Example 

The  theory  developed  will  now  be  illustrated  with  the  numerical  example.  The  algorithms  described  above 
have  been  implemented  by  the  computer  program  developed  in  Fortran  and  Visual  Basic  languages. 

The  blank  data  and  gear  machine-tool  settings  have  been  adapted  from  the  Gleason  summary  and 
represented  in  Tables  2  and  3. 

The  final  machine-tool  settings  for  pinion  generation  have  been  determined  as  the  output  of  the  computer 
programs  developed  and  are  shown  in  Table  4.  The  results  of  TCA  obtained  by  the  developed  computer 
programs  are  shown  by  Figs.  23,  24,  25,  26,  and  27  for  aligned  and  misaligned  gear  drives. 

The  investigation  accomplished  confirms  that  the  longitudinal  path  of  bearing  contact  developed  is 
stabilized  and  the  predesigned  parabolic  function  absorbs  the  transmission  errors  caused  by  misalignment. 
The  magnitude  of  transmission  errors  caused  by  misalignment  might  be  controlled  and  does  not  exceed  7 
arc  seconds  in  the  example. 


Table  2:  Blank  Data 


Applied  parameters 

Pinion 

Gear 

Number  of  teeth  V, , 

10 

52 

Shaft  angle  y 

90° 

Mean  spiral  angle  /?j ,  /J2 

32°30' 

32°30’ 

Hand  of  spiral 

Right  hand 

Left  hand 

Whole  depth  (mm) 

16.637 

16.637 

Pitch  angles  71,72 

10°53' 

797. 

Root  angles  7;.,  »7;.^ 

942' 

75°8' 

Face  angles  7/^,7/^ 

14°52’ 

80°18’ 

Addendum  (mm) 

11.811 

3.429 

Dedendum  (mm) 

4.826 

13.208 

Face  width  (mm) 

71.120 

Diametral  pitch  (1/inch) 

3.156 

Clearance  CpC2  (mm) 

1.400 

Mean  cone  distance  (nun) 

177.521 
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Table  3:  Gear  Machine  Tool  Settings  and  Installment  of  Head-Cutter 


Applied  parameters 

Data 

1 .  Blade  Specifications 

Blade  angle  at  mean  point  M 

20” 

Parabolic  coefficient 

0.002 

Cutter  Radial  R,,  (mm) 

177.800 

Cutter  point  radius  (mm) 

177.800  +  2.4123 

Blade  edge  radius  p„,  (mm) 

2.794 

Point  width  (mm) 

4.826 

Gear  Machine-Tool  Settings 

Machine  root  angle 

7610 

Machine  center  to  back  (mm) 

-2.007 

Horizontal  if  2  (™tt) 

100.152 

Vertical  (mm) 

137.897 

Table  4:  Pinion  Machine-Tool  Settings  and  Blade  Data 


Applied  parameters 

Concave  Side 

Convex  Side 

Blade  angle 

17 

23° 

Cutter  point  radius  i?^  (mm) 

174.496 

185.151 

Radial  setting  5^  (mm) 

150.301 

195.018 

Installment  angle 

62”3'45" 

61°32i2" 

Machine  center  to  back  (mm) 

-4.530 

8.134 

Sliding  base  (mm) 

0.339 

-1.795 

Blank  offset  (mm) 

19.947 

-20.314 

Ratio  of  cutting 

4.853 

5.773 

Machine  root  angle  Ym, 

9.7° 

ON 

Coefficient  of  modified  roll 

0.0302 

-0.0527 

Coefficient  of  modified  roll 

-0.0599 

0.1157 
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9.  Conclusion 


Based  on  the  research  performed  the  following  conclusions  can  be  drawn: 

(1)  The  developed  approach  of  design  has  been  successfully  applied  and  enabled  to  design  formate  cut 
spiral  bevel  gears  with  a  stabilized  bearing  contact  and  favorable  shape  of  transmission  errors  of  low 
magnitude. 

(2)  The  theoretical  results  obtained  have  been  confirmed  by  simulation  of  meshing. 
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Appendix  1: — Relationships  Between  Principal  Curvatures  and 
Directions  of  Mating  Surfaces 


Introduction 

The  procedure  of  local  synthesis  requires  the  knowledge  of  principal  curvatures  and  directions  of 
contacting  surfaces.  In  the  case  of  design  discussed  in  this  paper,  the  pinion  tooth  surface  is  represented  by 
three  related  parameters  and  the  determination  of  pinion  principal  curvatures  and  directions  is  a  complex 
problem.  The  solution  to  the  problem  is  simplified  by  representation  of  pinion  principal  curvatures  and 
directions  in  terms  of  principal  curvatures  and  directions  of  the  head-cutter  and  parameters  of  motions  of 
the  generating  process  (proposed  in  [8,  13]). 

Important  relations  between  principal  curvatures  and  directions  of  mating  surfaces  in  point  contact 
(proposed  in  [8, 13])  are  applied  for  local  synthesis  (see  section  6). 

Henceforth  two  types  of  instantaneous  contact  of  meshing  surfaces  will  be  considered:  (i)  along  a  Une,  and 
(ii)  at  a  point.  Line  contact  is  provided  in  meshing  of  the  surface  being  generated  with  the  tool  surface. 
Point  contact  is  provided  for  the  generated  pinion  and  gear  tooth  surfaces.  The  determination  of  the 
required  relationships  is  based  on  the  approach  proposed  in  [8,  14],  The  basic  equations  in  the  approach 
developed  are  as  follows 


(2)  (1)  (12) 

V,  =v,  -l-v 

(72) 

.  (2)  .  (1)  (12) 
n^  =n^  -Ko  xn 

(73) 

(74) 

dr  ■' 

Equations  (72)  and  (73)  relate  the  velocities  of  the  contact  point  and  the  tip  of  the  unit  normal  in  their 
motions  over  the  contacting  surfaces.  Equation  (74)  represents  the  differentiated  equation  of  meshing. 
Equations  (72)-(74)  yield  a  skew-symmetric  system  of  three  linear  equations  in  two  unknowns  x,  and 
of  the  following  structure 


a^x^+alX^=al^  (i=l-3)  (75) 

Here:  x,  and  Xj  are  the  projections  of  the  velocity  of  the  contact  point  in  the  motion  over  one  of  the 
surface  on  the  principal  directions  of  the  mating  surface.  In  case  of  line  contact  of  surfaces  the  solution  for 
the  unknowns  is  indefinite  and  the  rank  of  the  system  matrix  of  the  linear  equations  is  one.  In  case  of  point 
contact  of  surfaces  the  solution  for  the  unknowns  is  definite,  and  the  rank  of  the  system  matrix  is  two.  The 
properties  mentioned  above  are  used  for  the  derivation  of  the  sought  for  relationships  between  the  principal 
curvatures  and  directions  of  the  meshing  surfaces. 

Meshing  of  Surfaces  Z,  and  Zj :  Surfaces  Z,  and  Zj  are  in  point  contact  and  their  meshing  is  considered 
in  fixed  coordinate  system  S,  (Fig.  20).  Equations  (72)-(74)  yield  the  following  system  of  three  linear 
equations  [8, 14] 


(1)  ,  (1)  ..  , 


(76) 


where 
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(77) 


(II  (I) 

V  =  V  •  e 

*  s  r  s 


•e 


Q 


,  ,  2  (12)  ,  .  2  (12) 
a^^=k^-kj- COS  a  -kfj  sin  a 

(12) 

^12  =^21  =0.5(/:^  - kfj) sin2<7 


,  (^2)  (12)^  ^ 

^13  =^31  =“^'5  ^'5  +[(nxw  )*ej 

,  ,  .  2  (12)  ,  2  (12) 

^22  =k^-kj-  sin  <7  “  K;,  cos  <7 

(12)  (12) 

^23  =  <^32  =  “V'9  +[(nX<0  )•«,] 

,  ,  (12K2  ,  ^  (12)^2  (12)  (12) 

^i3  =  k,{v,  )  +fc/y^  )  -[(nxo  )-v  ] 

-n-[((t)"*xv|^*)-((o'^*xv|i’)]+m2,(nxk2)-r 


(78) 


The  following  are  considered  as  being  known:  point  M  of  tangency  of  surfaces  Z,  and  Z2 ,  the  common 
surface  unit  normal,  the  relative  velocity  ,  the  principal  curvatures  and  and  directions  and 
Cq  on  Z2  at  M,  and  the  elastic  deformation  5  of  surfaces  at  M.  The  goal  is  to  determine  the  principal 

curvatures  and  k,,  and  the  angle  formed  by  the  unit  vectors  e,  and  e, . 


The  velocity  v*/*(/=7,  2)  of  the  contact  point  on  surface  Z,  has  a  definite  direction  and  therefore  the  rank 
of  the  system  matrix  (76)  is  one.  This  property  yields  the  following  relation: 


'll 

fll2 

«13 

'21 

^22 

^23 

'31 

dyi 

33 

(12)  , 

F{kf,k^,k^,k  ,a  ,/n2|)  =  0 


(79) 


The  sought  for  solution  for  kj  ,  k,,  and  (T^'^*  can  be  obtained  if  the  following  parameters  will  be  chosen: 

the  derivative  ;  the  ratio  a!  5  ,  where  a  is  the  major  axis  of  the  contact  ellipse;  the  direction  at  M  of  the 
tangent  to  the  contact  path  on  one  of  the  contacting  surfaces  Z,  and  Z2.  The  relation  between  the 
directions  at  M  of  the  tangents  to  the  contact  paths  on  both  surfaces  is  represented  by  the  equation  [8, 14] 

.  „  -"13^?  +(«33  +  ai3''s  )tanT72  ,oo^ 

tan7j,= - 2 -  - - - -  (80) 

^33  ^23 tQXlTl2) 


Choosing  r]^  at  point  M,  we  can  determine  r|^ . 

(12) 

Procedure  of  Determination  of  kj^ ,  kf^  and  cr 

Step  1:  Determine  r|^  choosing  r]2  • 

Step  2: 

(1)  a33 

V,  = - - 

aj3  +  a23tanTJi 


(81) 
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(82) 


Step  3: 


Step  4: 


Step  5: 


Step  6: 


where  82=^-^- 
Step  7: 


Step  8: 

where  • 

Step  9: 

Step  10: 


(1)  _  a^^t^^nr|^ 

“ 


^13+^23 


_5 

2 

a 


A  =  - 


K. 


‘23 


4A' 


/  (‘^2  ^  .  (1)  2 

(V,  )  +(v,  ) 

(1)  ,  (1) 

+^23V,  ,  24 

,  (1),2  (1)  2 

(v,  )  +(v,  ) 


^22 


(1).2  (1)  2 

(n  )  +(v,  ) 


(1) 

S 

(1) 


U)^.  0)^2 

■«23^  +(V,  ) 

(1)  (1)  (1)^ 

9  2- 


“■^13^5 


0)  ,  (1)^2 

+  +(V^  ) 


^23^  q 


tan  2(7 


^^12 _ 

82  “(^11  ^^22) 


^i=- 


sin  2(7 


(12) 


(1)  _  ^(2)  _ 

Aj;  —  Aj-  Aj- 


h={K'^'-g,)l2 


(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 

(90) 


The  procedure  provided  above  can  be  used  to  obtain  the  sought  for  principal  curvatures  and  at  point 
M  of  tangency  of  surfaces  E,  and  Ej  and  the  principal  directions  on  E,  at  M. 

Meshing  of  Surfaces  E,  and  Ep ;  The  tool  surface  Ep  generates  the  pinion  tooth  surface  Ej .  Surfaces 
El  and  Ep  are  in  line  contact  and  point  B  is  the  given  point  of  the  instantaneous  line  of  contact.  The 
meshing  of  surfaces  is  considered  in  S„,  .  At  point  B  the  following  are  assumed  as  given;  the  curvatures  k^ 
and  k^  of  siuface  Ep  ;  the  unit  vectors  and  of  principal  directions  on  E, ;  the  surfaces  unit  normal ; 

the  relative  velocity  v***' .  The  goal  is  to  determine  the  principal  curvatures  k^,  and  k,  of  surface  Ep ,  and 
the  angle  cr^’^’that  is  formed  by  the  unit  vectors  and  Cp  . 


Equation  (72)-(74)  yield  a  system  of  three  linear  equations 
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(91) 


<-,vr+4vr’=<-: 


(pi 


(i=l  -3 ) 


where 


(p)  (p)  (p)  (p) 

V/  =v,  •e^,v,  =.v,  -e,, 


(92) 


The  direction  of  is  indefinite  since  surfaces  and  Zj  are  in  line  contact.  Therefore,  the  rank  of 
system  matrix  of  equations  is  equal  to  one.  Using  this  property,  the  following  equations  are  obtained: 


Here: 


where 


tan2<T^‘^U- 


-2d\zd2i 


‘^23“^13~(^/  ~^/j)^33 


0  2 
<(f3  +  ^23 

<^33 


7  .  OPi  ('P\  1 

Jjj  =  +  [(n  X  (j)  )  •  ] 

7  7  (12)  r/  "P'x  n 

^23  =  “^*'^*  +[(nxc)  )-eJ 

)  “^*(‘'/>  )  -[(nxo)  )-v  ] 

r/  Cl  <P>x  /  <P>  (‘\l 

-n  [((0  xv„  )-(co  XV,,)] 


(1/j)  (Ip)  (Ip)  (Ip) 

Vf  •«/  .  V,  =V  -e,, 


(93) 


(94) 


(95) 
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Vh 

Figure  1. — ^Illustration  of  tangency  of  surfaces  Xj  and  X2- 
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(  b  )  Convex  side 


(  c  )  Concave  side 


A, 


Figure  6. — Illustration  of  gear  cutter  blades  and  generating  surfaces. 
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Convex  side 


Concave  side 
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Figure  11. — Geometric  model  of  pinion  tooth  surfaces. 


Figure  12. — ^Illustration  of  parameters  rj2  and  a  applied  for  local  synthesis. 
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Local  synthesis 


Pinion  machine-tool  settings 


5  Pinion  tooth  surface  Gear  tooth  surface 


Simulation  of  meshing  and  contact 


change 


Straight  line  projection 
Lt  of  path  of  contact?^ 


8  Output  pinion  machine-tool  settings 


Stop 


Figure  14. — Flow  chart  for  procedure  1:  design  of  bearing  contact. 
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Figure  15. — ^Illustration  of  computations  for  determination  of 


C  End  ) 


Figure  16. — Row  chart  for  procedure  2, 
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Figure  18. — ^Flow  chart  for  procedure  3. 
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Figure  21. — For  derivation  of  pinion  machine  tool  settings. 
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Figure  23. — TCA  results  output  by  Visual  Basic:  simulation  of  meshing  without  misalignment. 
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Figure  24. — ^TCA  results  output  by  Visual  Basic:  simulation  of  meshing  with  misalignment  AH  =  0.1  mm. 


NASA/CR— 2001-210894 


50 


Figure  25. — TCA  results  output  by  Visual  Basic:  simulation  of  meshing  with  misalignment  AD  =  0.1  mm. 
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Figure  26. — TCA  results  output  by  Visual  Basic:  simulation  of  meshing  with  misalignment  AE  =  0.1  mm. 
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Figure  27. — TCA  results  output  by  Visual  Basic:  simulation  of  meshing  with  misalignment  Ay=  3'. 
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